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Abstract 

In this article the Butler-Volmer equation used in describing Lithium- 
ion (Li-ion) batteries is discussed. First, a complete mathematical 
model based on a macro-homogeneous approach developed by Neuman 
is presented. Two common mistakes found in the literature regarding 
a sign in a boundary conditions and the use of the transfer coefficient 
are mentioned. The paper focuses on the form of the Butler-Volmer 
equation in the model. It is shown how practical problems can be 
avoided by taking care in the form used, particularly to avoid difficul¬ 
ties when the solid particle in the electrodes approach a fully charged 
or discharged state or the electrolyte gets depleted. This shows that 
the open circuit voltage and the exchange current density must depend 
on the lithium concentration in both the solid and the electrolyte in 
a particular way at the extremes of the concentration ranges. 
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1 Introduction 


Lithium-ion (Li-ion) batteries have become very popular in the last years as 
efficient energy storage devices. A mathematical model demonstrating the 
role of the key factors in the battery operation can be very helpful for the 
design and optimization of new models and also for the real time control of 
its performance. 

Based on a macro-homogeneous approach developed by Neuman [1], sev¬ 
eral mathematical models have been developed for these purposes (cf. [2j[3j[3j 
El Eld El El HD], which include the main physics present in charge/discharge 
processes. A complete mathematical model is presented in Section [2] includ¬ 
ing a system of boundary value problems for the conservation of Lithium 
and conservation of charge in the solid and electrolyte phases. In that sec¬ 
tion we remark that several authors have included a boundary condition that 
is not correct (see m and Remark El) • Although the authors have probably 
used the correct boundary condition when solving the model numerically the 
reader should be aware of the right choice. A second common mistake in¬ 
volving the role of the transfer coefficient in the equation for conservation of 
charge in the electrolyte is also pointed out (see Remark |2j). 

Once we have presented the model, Section [3] is devoted to the Butler- 
Volrner equation. It is shown that care must be taken in the form of this 
equation especially in the limiting cases where the electrodes approach a fully 
charged or discharged state or the local electrolyte concentration reduces sub¬ 
stantially. In such cases it is necessary to ensure the equation allows particles 
to discharge or charge but only in direction that ensures concentrations re¬ 
main in the physically relevant regime. 


2 Mathematical model 

2.1 Generalities 

Let us assume we have a general binary electrolyte (i.e. a single salt composed 
of one kind of cation and one kind of anion) conducting electricity and there 
is no convection involved in the process. In order to give a general idea of 
the equations involved in the model detailed in Section [2721 below, we use the 
infinite dilute solution theory (we could also use the concentrated solution 
theory, see PQ, which is more complex and accurate but leads to equations 
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of the same type). 

The accumulation of ions is due to the net flux and the production, i.e. 

f)r- 

= —V • iVj + Ri (ie {+,-}), (1) 

where Ri (mole m -3 s _1 ) is the reaction rate and A) is the flux (mole m -2 
s^ 1 ) of anions (AL) and cations (N+) in the electrolyte. This flux is due to 
migration and diffusion and is given by 

Ni = -ZiUiFciV^e - AVq (i G {+, -}), (2) 

where c + ,c_ are the molar concentration (mol m -3 ) of cations and anions, 
z + , z_ are the number of protons charges carried by a cation and by an anion, 
u + ,u_ are the mobility of cations and anions (m 2 mole A 1 s -1 ) and D + , D_ 
are the diffusion coefficient of cations and anions (m 2 s” 1 ). 

We assume the solution is electrically neutral, i.e. 

z + c + + z_c_ = 0 (3) 

Let u +1 v_ be the number of cations and anions produced by the dissoci¬ 
ation of one molecule of electrolyte, i.e. 


z + is + + z_v_ = 0. 

and 

c e = — = — (so called, concentration of the electrolyte). 

is + u_ 

Assuming u + and il_ are constant, if we multiply dH) by A- we get 

8r /? 

- - z + u + FV • (cVA) - V • (D+Vc) = -±, 

Ot ZA|_ 

dc R 

— - z_u_FX7 ■ (cV0 e ) - V ■ (D_Vc) = —. 

ot V _ 

Subtraction of these two equaitons gives 

(. z+u + - z_u_)FV ■ (cV0 e ) + V • ((£>+ - D_)Vc) = --?*+, 

Z/_ u + 
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which can be used to eliminate V0 e from either of the previous equations in 
order to get the equation for the conservation of ions 

- V • (DeVCe) = R, (4) 

where 

z + u+D~ - Z-U-D+ Z + U +Tfi ~ Z - U -Tfi 

D e = -, and K = -—. 

Z + U + — Z-U- ' Z + U + — Z-U- 

The diffusion D e and reaction R coefficients represent a compromise between 
the diffusion and the reaction coefficients of the anion and the cation. 
Regarding the electrical current in the electrolyte, we have that 


i e = F(z + N + + z_N_), 

where i e is current density (A m -2 ). Then, from (fT j) - (11]) it can be deduced 
that 

V • i e R- R + , ^ 

-F =- ( 5 ) 

Z+P+f v_ p + 

and, using the Nernst-Einstein equation, D t = RTui , that 

_ , RT n ( f_|_ t—\ . , 

i e = -/cV0 e - — — H-V In c e , (6) 

F \z + z-J 

where n = F 2 (( z + ) 2 u + c+ + (z-) 2 U-C-) is electrical conductivity of the elec¬ 
trolyte and 


t, = 


z+u+ 


t- = 


-Z_U- 


= 1 -t: 


z + u + — Z-U- ' 


z + u + — Z-U- 


are the so-called transference numbers. Hence, from (J5J) and (J6|) we obtain 
the equation for the conservation of charge 


- V • (/cV0 e ) 


RT 

~Y 


— + —) V-(/evince) 

z + z-J 


z + v + F 




• (7) 


Remark 1 For Li + batteries with z+ — 1, z_ — — 1, v + — p_ — 1 and 
R- = 0 (i.e. the cation is the only ion reacting at the electrode) the following 
equations are satisfied in the electrolyte: 


dc e 

dt 


V • (-D e Vc e ) 


1-4 

F 


•Li 


(conservation of ions) 
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-V • («V0e) 



1) V • (kV In c e ) = j Ll (conservation of charge). 


We have used that, due to the Faraday’s laws of electrolysis, R + = yj h \ 
where j Ll = Qjj is the reaction current (A m~ 3 ) resulting in production or 
consumption of Li + . 


2.2 The complete model 

A typical Li-ion battery cell has three regions: A porous negative electrode, 
a porous positive electrode and an electron-blocking separator. In all these 
regions there is an electrolyte containing various charge species, including 
lithium, that can move all through the cell in response to an electrochemical 
potential gradient. 

A ID electrochemical model is considered for the evolution of the Li 
concentration c e (x , t ) (mol m -3 ) and the electric potential 4> e (x, t ) (V) in the 
electrolyte and the electric potential 4> s (x, t ) (V) in the solid electrodes, along 
the ^-direction, with x G (0, L ) and L = L 1 + 6 + L 2 being the cell width (m). 
We assume that (0,Li) corresponds to the negative electrode, (Li,Li + <f) 
corresponds to the separator and (Li + <f, L ) corresponds to the positive 
electrode. The lithium behavior in the electrolyte is coupled with a diffusion 
model for the evolution of the Li concentration c s (x',r,t ) in a generic solid 
spherical electrode particle. Spherical symmetry of this diffusion is assumed 
along the radial r-direction in the region r G [0, R.f\ where R s (m) is the 
average radius of a generic particle. Using a ID approximation for the entire 
battery is valid since the characteristic length scale of a typical Li-ion cell 
along the x-axis is on the order of 100 pm, whereas the characteristic length 
scale for the remaining two axes is on the order of 100,000 fim or more (cf. 
i). Note that R s can be different in each electrode {Rs~, R s ,+)• 

Based on the conservation equations deduced in Remark [T], and the mod¬ 
els appearing in the literature (cf. [21 3] @1 [5,12, T\ Ef [9, TO], along with the 
additional assumption of constant diffusion and activity electrolyte coeffi¬ 
cients, a model for the performance of a battery at constant temperature, is 
given by system of equations (1ST)— (1TTT) : 


5 


£e dt De dx ( £ ° dx) F + jLl ’ m ( °’ L) X ( °’ tend) ’ 
c)c 3c 

^ (M) = ^ (Lit) = 0 ’ 

C e (x, 0) = Ce i0 (x), xe(0,L), 


For each x G (0, L\) U {L\ + 5, L) : 

^ ( r ^) = °> 111 (°’^) X (°^end), 

^(x;0,t) = 0, * e (Mend), 

c s (x; r, 0) = c S)0 (x; r), 


For each t e (0. i en d) : 

( e ° k W) + ii ~ 2 *° in (ce) )= jU in (o ' l) - 

^(0, J ) = ^(L,t) = 0, 

ax ax 


For each t e (0,f e nd) : 


-e s a-^- = - j Ll in (0, Li) U (Li + 5,L), 
s,(0)a(0)^(0, t) = E,(L)a(L)^(L, t) = <U> 

Ms,. + \ _ Ms,. , A _ n 

1. ~g^ {Ll ' t) - ~a^ l ' Ll+li ' t) - 0 ' 

In this system of equations the independent variable are position x (m) and 
time t (s) with dependent variables c e = c e (x,t ) (mol m -3 ), c s = c s (x;r, t) 
(mol m -3 ), (p e = (j> e (x,t) (V), 0 S = cj> s (x,t) (V), and I = I(t) (A) the applied 
current. Parameters in the model are p the Bruggeman porosity exponent 
(nondimensional constant), D e the electrolyte diffusion coefficient (m 2 s^ 1 ), 
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t Q + the transference number of Li + , k = n (c e (x, t)) the electrolyte phase ionic 
conductivity (S m _1 ), and A (m 2 ) is the cross-sectional area (also the current 
collector area), There are parameters that take different values in different 
regions and these are £ e (e e ,- if x G (0,Li), £ e ,se P if x G (Li,Li + 5) and 
£ e ,+ if x G ( L\ + S,L)) which is the volume fraction of the electrolyte, D s 
(D s _ if x G (0,1a), D s + if x G (Li + S,L)) which is the solid phase Li 
diffusion coefficient (m 2 s -1 ), £ s (e s _ if x G (0 ,Li), s S:+ if x G {L\ + 5,L)) 
which is the volume fraction of the active materials in the electrodes, a s (cr S) _ 
if x G (0, Li), cr Si+ if x G (L\ + 5,L )) which is the electrical conductivity of 
solid active materials in an electrode (S m -1 ), and j Ll (A m -3 ) is the reaction 
current resulting from intercallation of Li into solid electrode particles. 

For j Ll the following Butler-Volrner equation is commonly used (cf. p, UJ 

El 0 El El HD]) 


3£ s (x) . 


•Li 


RJx) 


t 0 exp 


a a ,F 


~V 


exp 


-Ot c F 


3 = 


RT'J 1 \ RT 
if x G (0, Li) U (Li + 8,L), 

if x G (Li, Li + 5) 


( 12 ) 


(here, for the sake of simplicity, we have considered the solid/electrolyte 
interfacial film resistance to be zero and therefore is not included in the 
above equation), where T (K) is the temperature of the cell, a s (x) = 

(m -1 ) is the specific interfacial area of electrodes, to = io(x, c s , c e ) (A m -2 ) is 
the exchange current density of an electrode reaction , a a , a c (dimensionless) 
are anodic and cathodic transfer coefficients for an electrode reaction, and rj 
(V) is the surface overpotential (V) of an electrode reaction. 

Most of the parameters in the Butler-Volrner equation are agreed upon 
in models but there is considerable variation in the behavior of to and r). 
Typically 


V = 


4> s 

o 


U(x,Cs), if x G (0, Li) U (Li + 8 , L), 
if x G (L\, L\ + 5), 


where U(x,c s ) is the equilibrium potential (V) at the solid/electrolyte inter¬ 
face (i.e. the open circuit voltage, OCV, of each electrode). An empirical 
function is usually exploited for U taking different forms in the two electrodes 
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and dependent only on the local surface concentration of lithium with 


if x G (0, Li), 
if x G (Li + 5, L), 

The functions U-, U + are typically obtained from fitting experimental data 
and the constants c S)+;max , c Sj _ )max (mol m -3 ) taken as the maximum possible 
concentration in the solid positive and negative electrode, respectively. The 
exchange current io is taken in different forms including being constant and 
being 

• f k- Cg a (eg _, max - c s )“ a c“ c if x G (0, Li), , , 

\ k+ Cg a (c S)+imax - c s )" a c“ c if x G (Li + S,L), 

where k_,k + are kinetic rate constants (A m - 2 + 6 "a+3a c mol~ 2aa_ “ c ). In the 
final section we concentrate on the appropriateness of the functional form of 
U and io. 

Remark 2 The second term on the left hand side of system KTU) is often 
written in the literature, using (1 —(see /If, /3J/J or 2(1 — (see (ff, fflf, 

0 . 0 , a. mn, instead of (1 — 2), which is the term deduced in Remark 

m). 



Remark 3 This system of equations does not have uniqueness of solution 
(if (f> s (x,t) and f> e (x,t) are solutions then f> s (x,t) + c(f) and f> e (x,t) + c(t ) 
are also solutions, for any function c(t)). A way of avoiding that is to set a 
reference value of <j> s (x,t) or f> e (x,t) at some point x. For instance we can 
impose f> s (0,t) = 0 for any t G [0, t eri f\. Some results regarding the existence 
and uniqueness of solution can be seen in m 


Remark 4 After solving the above model, we can estimate the state of the 
charge of the negative electrode SOC_(t) and of the positive one SOC + (t), 
and the cell voltage V(t), at time t, by computing 


pRs. 


SOC.it) = 


U{R S -Y Jo 


2 c s(x] r ,t) 
r -drcLr, 

^s,-,max 


SOC+{t) 


(L «- Li — 5) (Rg-\ 


f L / 

' L\+8 J 0 


Rs ’+ 2 c s(x'i r,t) 

r -dr day 
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V{t) = 4> S (L, t) - 0 s (o ,t) - ^l(t), 
where there is a constant film resistance of Rf. 


Remark 5 In ffif, f6f. fTjj, j2J/, JUIf, authors use the following (incorrect) 
boundary conditions at x = 0 and x = L, instead of those presented in m 


, „ _ c 8<t>. _ I(t ) 

S ’ dx ( ’ } “ S ’ + + dx ( ’ } “ A 


(14) 


We note m passing that if this incorrect condition is used then it can be 
proved (see m) that the corresponding system of boundary value problems 
does not have any solution unless I(t ) = 0. 


3 The Butler-Volmer equation 

The Butler-Volmer equation in the mathematical model presented here has 
a general functional form widely used in the literature (cf. |2j 3j HJ (i. 7j, SJ 
El HD]), but it is necessary to take care in the detailed functions used for U 
and ij as explained below. We note that both io and U may vary with c s 
and c e throughout the relevant range of these concentrations but we are only 
interested in ensuring that the behavior is appropriate when the electrodes 
get close to being fully intercollated c s = c s ±. niax . are fully depleted c s = 0, or 
where the electrolyte is completely depleted c e = 0. To discuss the possible 
behavior we consider the case c s —» 0 and then indicate how this can be 
extended to the other cases. 

Consider therefore the case c s —» 0 while U tending to a constant and 
i 0 is given by (1T3|) so that i 0 — > 0. It might be expected that the Li flux 
out of the particle will therefore cease as c s —> 0 thereby preventing negative 
concentrations in the solid. However, this condition also implies that if we 
start with a depleted particle then no flux can ever enter the particle. Such 
a physically unrealistic situation should not be allowed in the model. 

Alternative formulations of ?' 0 and U can avoid this problem in various 
ways and it is always possible to simply chose these functions and then impose 
some switching logic to turn the flux on or off as required to avoid such 
problems. However, it may be preferable to have the Butler-Volmer condition 
designed, through the modeling of suitable physical mechanisms, to ensure 
such physically irrelevant situations cannot occur 
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Hence we would like to have the property that, as c s —» 0 the flux can only 
take negative values. Such a property can readily be achieved by ensuring 
that U —> oo as c s —> 0 as this will ensure that the negative exponential 
dominates the Butler-Volmcr condition. A common physical condition is to 
consider that, as c s —» 0 the forward (cathodic) should be of zeroth order in 
c s while the reverse (anodic) reaction should be first order. Such behavior 
can be readily achieved by taking i 0 and U to have the local form, when 
c s -> 0, 

/ n c \ — RT 

i 0 ps (k Cg A and U pa —-- In (k c s ), 

F(a a + a c ) 

where k is a positive constant (which can be different for each electrode and 
also can change inside non-homogeneous electrodes). Such a formulation near 
this extreme of the concentration will automatically ensure the flux cannot 
become positive at this limit of the concentration reduces, thereby avoiding 
negative concentrations, and the flux can be finite and positive in this limit 
so the particle can be charged from a completely depleted state. Note that 
a reverse reaction of order greater than one can be analyzed but this is not 
usually considered. 

Such local behavior of the Butler-Volmcr condition is included in some 
formulation but these need to include all possibilities including extremely 
low electrolyte concentrations. 

Therefore, a reasonable Butler-Volmer formulation would be one where 

i 0 (x, Cs, C e ) = i a (x, C s ) “a+«c i c (x, Ce, Cs) “a+«c (15) 

and 

_ RT f i c (x,c e ,c s ) \ 

F(a a + a c ) V *a(®, c s ) ) ' 

The functions i a and i c must be strictly positive except at the extremes of 
the concentrations and locally these functions must take the form: 

i) as c s —> 0 

*a(®, C s ) ~ K{x) C s , 

with kg.(x) = /c a ,+ or k a -, if x e (0, L\) or x G (Li + S, L), respectively (again, 
non-homogeneous electrodes could also be considered). 

ii) as c s y c Sj= i- jmax 



z c \x 5 c e5 c< 


k 


c,s \Xy C e j (C s? -|- jma x C s 
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and 

iii) as c e —* 0 

i c (x,c e ,c s ) « k c<e (x, c s ) c e . 

Note there may also be a need to consider a maximum electrolyte concen¬ 
tration to avoid precipitation and this might be accommodated in a similar 
manner. 

Notice four important features of the approach proposed here in contrast 
with what is commonly found in the literature: 


1. The function U, the OCV, also depends on c e (instead of being inde¬ 
pendent of it). 

2. If c s —» 0, then i 0 —>• 0 and U —* oo (instead of tending to a finite value) 
so that 

i Ll ~ (W s exp s - </>e)^ - -Di exp ^-^^(0 S - 0 e )^ 

for some strictly positive constants Hi an D\. 

3. If c e —* 0, then i 0 — > 0 and U —> —oo so that 

jLl ~ 7T ( Z?2 6XP ^ “ D 2Ceexp - 0 e )^ 

for some strictly positive constants B 2 an D 2 . 

4. If c s —» c S) ± )inax , then z 0 —> 0 and U —> —oo and 

jLi “lK B3exp (^ w>, “ <w ) 


-D 3 (c s - c s , ±imax ) exp 

for some strictly positive constants B : > an D 3 . 


—oc c F 
R T 


(0s 0e 


These conditions ensure that the lithium flux is constrained to prevent con¬ 
centrations in the solid being taken into an unphysical regime and also ensur¬ 
ing that, if the solid is near such extreme conditions, the flux can move the 
system away from these. This ensures the battery can readily charge from a 
completely depleted state or discharge from a fully charged state. 
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Remark 6 The most simple case of functions i 0 and U satisfying the above 
conditions (of course, many other possibilities can be also considered) is ob¬ 
tained with 

■ i \ _ f K- C s if X E (0, Lx), 

* a(x,Cs) i K+c s ifxe(Lx + 5,L), 

■ / \ _ / ^c,-^e(^s,-,max Cs) if X (z (0, -^l), 

Mx,c s ,c eJ -| kc+Ce{Cs + ^ _ Cg) ifxe (L 1 + 5,L). 

Therefore, according to nT\) . 


U(x, c s , c e ) 


RT 

F(a a +a c ) 

RT 

F(a a +a c ) 


In 

In 


^c, — c e(Cs, —,n 


^a,— Cs 
^c, + C e (c s ^ )m 


^a, + C B 



if x e (0 ,Lx), 
if x e ( Lx + 6, L). 


Furthermore, according to / fldl) . 


io{x, c s , c e ) 


Q a 

k_Ce a+ac 


Q a 


k + c 


e 


(^s,-,max 

(^s,+,max 



Q=c 

n a a + a c 
^S 


C s )' 


c^a+c^c 
^S 


x e (o,Li), 
if x G (Lx + h, L), 


with k_ = kf c _kf & _ and k + = kf( + kff_, which coincides with KTR) when a a + 
a c = 1. 


Remark 7 Some authors (cf. mmm)) take a constant value for i 0 and 
some of them (cf. claim that it exhibits modest dependency on electrolyte 
and solid surface concentration. Although this can be valid for appropriate 
particular cases, in a general situation this does not seem to be valid, since 
i 0 may vary importantly, with extreme cases when the battery is either fully 
charged or fully discharged. 


4 Conclusions 

The Butler-Volmer equation is commonly used in the literature in order to 
take into account the electrochemical reactions that take place in battery 
electrodes. It is widely used in mathematical models for the simulation of a 
Lithium-ion battery, based on a macro-homogeneous approach developed by 
Neuman [1]. It has been shown that the way this equation is commonly used 
in these models may create situations where the system becomes non physical. 
Conditions on the functional form of the Butler-Volmer equation that ensure 
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this will not occur have been presented. Furthermore, some common mistakes 
that can be found in the literature regarding the mathematical equations of 
the Neuman-type models have also been pointed out. 
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